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U . ABSTRACT 

■ We present new spectroscopic and photometric data on the early-type overcontact binary 
TU Muscae. The analysis of the spectroscopic data shows that the line of sight to the system 

. . crosses three kinematically sharp and well-separated interstellar reddening sources and that the 

rS ' stars rotate synchronously. We present new radial velocities that are in good agreement with 

■ earlier optical velocities and, thus, do not confirm the systematically smaller velocities obtained 
from lUE spectra. The optical velocities are analyzed simultaneously with the photometric data 
to derive accurate absolute dimesions for the binary components. The results show that TU Mus 
consists of an 07.5 primary with Mi = 23.5±O.8M0, Ri = 7.48±0.08i?o and an 09.5 secondary 
with M2 — 15.3±0.4Mo, R2 — 6.15±0.07i?Q in an overcontact configuration and that the orbital 
period has remained constant over the three decades covered by the observations. These results 
might imply that the mass transfer seen in late-type overcontact binaries does not occur in their 
early-type counterparts. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: individual (TU Mus) — stars: 
main sequence 

1. Introduction Gr0nbech (1975, hereafter AG75) who list refer- 

ences to earlier work on this early-type overcontact 
The variability of TU Muscae (HD 100213) (Wilson 2001) binary. Given the paucity of eclips- 

was discovered by Oosterhoff (1928) and the sys- j^g, double-hned 0-type binaries, TU Mus has re- 

tem was first studied in detail by Andersen & ^^^^^^ surprisingly little attention from observers. 
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AG75 published fully covered uvby light curves 
and fourteen radial velocities for each star near 
quadratures from photographic spectra of moder- 
ate resolution. Stickland, ct al. (1995, hereafter 
S95) presented twenty-four radial velocities for the 
primary star and twenty-three for the secondary 
from high resolution International Ultraviolet Ex- 
plorer (lUE) spectra. TU Mus was bright enough 
(Vt = 8.4) to be observed by Hipparcos and a light 
curve with about a hundred points was published 
in the Hipparcos Epoch Photometry Annex (ESA 
1997). Two analyses of these data using modern 
light curve synthesis tools have been published: a 
preliminary report on the current work by Terrell 
(2002) and an analysis of the AG75 photometry 
by Wilson & Rafert (1981). 

The two sets of radial velocities give quite 
discrepant results for the absolute dimensions. 
AG75 obtained photographic spectra with a dis- 
persion of 20 A mm and found Ki=25il.3 ±4.8 
km sec-i and ii:2=370.5 ±4.3 km sec"^. S95 
found ifi=216.7 ±2.7 km scc-^ and A'2=345.4 
±3.1 km sec~^ from their lUE spectra that had 
signal-to-noise ratios (S/N) of 10-20. Adopting an 
inclination of 76° , these numbers result in masses 
of Ml = 17.2M0 and Ma = IO.8M0 for the lUE 
velocities and Mi = 23. SM© and M2 = 15. SM© 
for the optical velocities, a difference much larger 
than the error limits allow. In order to explore 
the nature of this discrepancy, we obtained high- 
resolution, high S/N spectra of TU Mus near the 
quadratures. Although we have a limited num- 
ber of spectra due to the difficulty of obtaining 
observing time on large telescopes for the study 
of eclipsing binaries, we find very good agreement 
with the AG75 results. 

Given the good agreement between our veloci- 
ties and those of AG75, we performed a simulta- 
neous solution of those velocities, the AG75 uvby 
photometry, the Hipparcos photometry and new 
BV photometry from the 2002 observing season. 
We used the 2003 version of the Wilson-Devinney 
(WD) program (Wilson & Devinney (1971), Wil- 
son (1979), Wilson (1990)) and present new ele- 
ments for the binary, showing that the period has 
remained remarkably constant over the roughly 
thirty year span of all the observations. 



2. Spectroscopy 

Nine ESQ 2.2m + FEROS spectra have been se- 
cured, in three groups of three spectra each (Table 
1). The wavelength range of the spectra extends 
from 3900 A to 9200 A with R = 48, 000 and a 
400 second exposure time. The S/N is evaluated 
around H/3 and the orbital phase is computed with 
the ephemeris given in section 4. 

2.1. Rotational velocity 

AG75 estimated the rotational velocities of the 
primary and secondary stars in TU Mus as 285 
and 240 km sec~^ (±10%), while S95 derived sig- 
nificantly slower rotations: 250 and 195 km sec~^ 
(±10%), respectively. To derive the rotational ve- 
locity we used the relation 

Vrot = 42.42 X HIW - 35km sec"^ (1) 

calibrated by Munari & Tomasella (1999) on the 
width at half maximum of He I 5876 A in high 
resolution spectra of O and B stars. The decovo- 
lution of the He I 5876 A line profile in the two 
TU Mus quadrature spectra gives a HIW of 7.26 
and 7.35 A for the primary, and 6.37 and 6.36 A 
for the secondary, which corresponds to 273 and 
277, and 234 and 236 km sec~^ respectively. The 
Munari and Tomasella relation was calibrated us- 
ing single, and thus axially symmetric, stars. The 
components of TU Mus depart quite strongly from 
axial symmetry so we investigated its influence on 
the HIW using the line profile feature of WD. We 
computed the HIW for TU Mus and for a spheri- 
cal star with radius equal to the r^acfe radii of the 
TU Mus components. We found that the HIW for 
TU Mus was about 5% smaller than for the spher- 
ical star. Applying this correction to our data, we 
adopt 290 and 248 km sec~^ as the rotational ve- 
locities of primary and secondary. Of course, the 
concept of a rotational velocity is of less obvious 
usefulness for stars like TU Mus than for axially 
symmetric ones. The question arises as to where 
on the equator of TU Mus the above rotational 
velocities apply. Since our spectra were taken at 
quadrature phases, rback seems appropriate and 
we find values of 291 and 242 km sec~^ there from 
our light curve solution (which assumed synchro- 
nism), in excellent agreement with the line profile 
results. Our data support the AG75 values and 
rule out the smaller values measured by S95. 
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Table 1 

ESO 2.2m + FEROS Spectra of TU Mus 



Spectrum Number 


Date 


UT of mid-exposure 


S/N 




Orbital Phase 


1346, 1347, 1348 
1629, 1630, 1631 
1693, 1694, 1695 


Feb 22, 2003 
Feb 25, 2003 
Feb 26, 2003 


07:08:56, 07:18:16, 07:27:31 
06:07:24, 06:16:39, 06:25:59 
00:45:27, 00:54:47, 01:04:07 


75, 68, 
160, 165, 
125, 130, 


72 
158 
120 


0.032 
0.211 
0.771 



To confirm the MT method we have appUed 
it to the He I 5876 A fine of the B1.5Vn star 
HIP 77635 that wc observed under identical in- 
strumental conditions soon after TU Mus on 
Feb. 22. Its rotational velocity is reported as 
306 km sec~^ in the Bernacca & Pcrinotto (1970) 
catalogue of rotational velocities. The HIW of the 
line turned out to be 8.16 A which corresponds 
to 311 km sec""'^ according to Equation 1. This 
excellent agreement supports our confidence in 
the rotational velocities we have measured for the 
components of TU Mus. We have also computed 
synthetic, broadened profiles for He I 4471, 6678 
and 7065 A and the agreement with the observed 
profiles is also excellent. 

2.2. Distance 

Assuming a luminosity ratio L2/-^i=0.55 from 
our light curve solution, the Vt=8.41 magni- 
tude transforms into Johnson V magnitudes of 
Fi=8.88 and ^2=9. 52, following Bcssell (2000). 
With i!^B_y=0.19 (see below) and assuming an 
07.5V classification for the primary, the distance 
to TU Mus amounts to 4.77 kpc using the absolute 
magnitudes calibrated by Houk, Branch & Win- 
frey (1994) (Afy=-5.1, -4.5 for 07.5V, 09.5V 
stars) and a ratio of total to selective absorption 
of Rv =3.1. The same calculation for the 09.5V 
secondary results in a distance of 4.86 kpc. We 
therefore adopt 4.8 kpc as the distance to TU 
Mus. 

2.3. Reddening 

The interstellar lines of Na I offer a means of 
estimating the reddening affecting TU Mus inde- 
pendent of the knowledge of the intrinsic colors. 
The profile of Na I Dl and D2 lines on our FEROS 
spectra of TU Mus are shown in Figure 1. They 
are clearly composed of three separate compo- 
nents. Their average position, width and equiva- 



lent width have been derived from multi-Gaussian 
fitting on all 9 available spectra, and the results 
are summarized in Table 2, where the errors are 
the standard errors. The 9 individual spectra gave 
extremely consistent results. The multi-Gaussian 
m,ean fit of Table 2 is overplotted on the D2 line of 
spectrum #1629 in Figure 1 to show the accuracy 
of the fit. 

The equivalent widths of the three components 
have been transformed in Table 2 into the corre- 
sponding color excess using the relation calibrated 
by Munari & Zwitter (1997), with a total color 
excess amounting to ii^B-y=0.19±0.01. Similar 
profiles are observed for K I 7698 A and H & K 
Ca II interstellar lines. The Na I widths in Table 
2 are not corrected for instrumental width, which 
amount to 5.0 km sec~^ from unresolved telluric 
lines in the vicinity of Na I lines (thus intrinsic 
width does not exceed 3 km sec~^). The line of 
sight to TU Mus therefore crosses three kinemati- 
cally very sharp and well separated sources of red- 
dening, of pure interstellar origin. In fact none of 
the components show radial velocity change with 
orbital phase and none share the -4 km sec~^ sys- 
temic velocity. 

Of interest is the heliocentric radial velocity of 
the three components. The accuracy reported in 
Table 2 is not fictitious because nearby telluric 
absorptions show the same wavelength stability 
from one spectrum to the other, and it matches 
the expected high FEROS spectrograph stability 
suitable for extra-solar planet searches. The lim- 
ited resolution of the AG75 spectra did not allow 
them to recognize the multi-component structure 
of the interstellar lines (H and K lines of Ca II), 
for which they estimated a — 7.3±0.8 km sec""*^ ra- 
dial velocity. Averaging the three components in 
Table 2 weighted according to equivalent width, a 
radial velocity of — 0.7±0.2 km sec^^ is derived for 
an unresolved interstellar Na I profile. However, 
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Table 2 

Components of the Interstellar Lines in the TU Mus Spectrum 



Quantity 


Component 1 


Component 2 


Component 3 


Total 


RV (km sec-i) 


-15.51±0.09 


+5.8±0.2 


+ 15.6±0.1 




Width (km sec"^)*" 


6.00±0.04 


5.8±0.1 


5.3±0.2 




Na I D2 E.W. (A) 


0.225±0.002 


0.216±0.002 


0.099±0.003 






0.080 


0.075 


0.031 


O.lQibO.Ol 



* Uncorrected for the dominating ~5.0 km sec ^ instrumental width 




58B9 5890 5895 



heliocentric wavelength (A) 

Fig. 1. — The interstellar lines of Na I D 

(spectrum #1629) with superimposed the three- 
component Gaussian fitting of Table 2. 

a —7.2 ± 0.4 km sec~^ value averaged by equiva- 
lent width is obtained on our spectra for both H 
& K Ca II as well as K I interstellar lines. While 
showing exactly the same two major components 
of Na I Dl & D2, they fail to show the third, red- 
dest one. 

The radial velocities of the three components 
are plotted in Figure 2 where they are compared 
with the expected values of the mean galactic ro- 
tation along the Hue of sight to TU Mus (^=284. 8, 
6=-4.1 deg) from the Hron (1987) formalism with 
A=-'17 km sec^^ and km sec~^ pc~^ and 

Pont, et al. (1994) values for {uo, Vo, Wo)- It 
is evident that the mean Galaxy rotation curve 
does not reproduce the expected velocities up to 
the ^4.8 kpc distance to TU Mus (note that the 
line of sight to the target is still mostly within the 
2: = ±150pc dust layer around the galactic plane 



at the estimated ~4.8 kpc distance). The dis- 
crete velocity field maps of Brand & Blitz (1993), 
based on H II regions and reflection nebulae, do 
not perform much better, which could be antici- 
pated given the scant and widely spaced data used 
in their analysis in that part of the Galaxy. Only 
the stronger component can be accounted for and 
located at ~2.5 kpc distance by the Brand & Blitz 
(1993) map. The limited individual reddening im- 
plied by the equivalent widths of the three inter- 
stellar lines suggests that only a few limited and 
kinematically scattered bubbles of the interstellar 
medium are crossed by the line of sight, with a 
diffuse component playing a marginal role, if any. 
The kinematic information from the interstellar 




distance from the Sun (kpc) 



Fig. 2. The radial velocity of the three com- 
ponents of Na I interstellar lines (arrows) com- 
pared with behaviour of the general galactic ro- 
tation curve of Hron (1987, continous line) and 
the local mapping by Brand and Blitz (1993, solid 
points) up to the distance of TU Mus (vertical 
dashed line). 
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lines cannot therefore be used to constrain the dis- 
tance to TU Mus. 



2.4. Radial velocities 

The radial velocity of TU Mus has been mea- 
sured on the two groups of spectra in Table 1 ob- 
tained close to quadrature phases. Those around 
orbital phase 0.03 are very close to conjunction 
and thus have insufficient velocity separation be- 
tween the two components to obtain reliable ra- 
dial velocities. Line splitting at quadrature is 
wide and allows an easy and firm determination 
of the radial velocity, as shown for H/3 in Figure 
3. The radial velocities have been obtained by 
Gaussian de-convolution of the profiles into two 
components, which provide an excellent overall fit 
to the observed line profile. The measured lines 
are He I 4016.218, 4471.507, 4921.929, 5015.675, 
5875.651, 6678.149, 7065.276 A; He II 4199.831, 
4541.589, 4685.682, 5411.524 A; and H7, H/3, Ha. 
The mean of the radial velocities (and associated 
standard error) derived from these 14 lines is in 
Table 3. 

Our high dispersion, high S/N radial velocities 
are very close to the photographic ones by AG75 
and definitively rule out the smaller ones obtained 
by S95 on low S/N lUE SWP high resolution spec- 
tra. This is clearly demonstrated by Figure 4 that 
shows the He I 5876 A line of spectrum ^ 1629 su- 
perimposed with the two component fit using our 
radial velocities in Table 3 and the corresponding 
values for AG75 and S95. 

2.5. Spectral classification and metallicity 

Both stars are 0-type, showing the whole series 
of expected He II lines. Their spectral types are 
somewhat different, as the comparison of H/3 and 
He I 4471 A profiles in Figures 3 and 4 suggests. 
The ratio of the two H/3 lines is 1.8 while those 
of the He I is 1.3, favoring a two spectral subtype 
higher temperature for the primary. A reasonable 
estimate seems to be 07.5 for the primary and 
09.5 for the secondary. 

Other than the hydrogen and helium lines, only 
two C III lines are clearly visible at 4650 A and 
4069 A(Si IV lines are unsuitable in this regard.) 
Both C III lines are actually blends of three indi- 
vidual lines with complicated profiles caused by 
wide wavelength separation within the triplets, 



1.1 W 





1. I .1'/ 



1860 -1B70 4880 

Lric WHvelengLl, (A) 



Fig. 3.- The H/3 profile of TU Mus for the aver- 
aged Feb 25 and Feb 26 spectra of Table 1. The 
ticks mark the primary and secondary radial ve- 
locity from Table 3. 




4465 4470 4475 44B0 

heliocentric wavelength (A) 

Fig. 4.— The Hel 4471 A profile for Feb 25 spectra 
of Table 1, superposed with the double Gaussian 
fitting with our velocities from Table 3 (thick solid 
line), corresponding AG 75 velocities (thin solid 
line) and S95 ones (dashed line). 
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Table 3 

Heliocentric radial velocities of TU Mus 



Spectrum Number Orbital Phase Primary Secondary 

This Paper AG75 S95 This Paper AG75 S95 



1629-30-31 0.211 -244 ±5 -246 -225 +350 ± 5 +364 +324 

1693-94-95 0.771 -1-237 ±4 -1-249 H-197 -380 ±8 -362 -348 



fast rotation, and the SB2 nature of the spec- 
tra. However, a comparison with synthetic spectra 
does not support a carbon content, and therefore 
metaUicity, far from the solar value. 

3. Photometry 

TU Mus was observed by RHN on 10 nights 
from 2002 April 13 to Jimc 4 at the Mount John 
University Observatory while he was a guest of 
the University of Canterbury in Christchurch, New 
Zealand. The f/13.5 Cassegrain focus on the 0.61 
meter Optical Craftsmen telescope was used to- 
gether with a Santa Barbara Instrument Group 
ST-9E CCD camera, utiUzing the KAF-0261E 
chip (512x512) with 20 micron pixels. A tele- 
compressor lens operated at 2.11x compression re- 
duced the focal ratio to 6.4 and gave a field of view 
of 9' X 9'. The B and V (Johnson) filters were by 
Schuler. 

Reduction of the CCD images and aperture 
photometry was done with the MIRA software 
package.^ Bias and dark removal phis flat field- 
ing (dome flats) were done in the usual way. 
The comparison star was GSC 8984:0579 and 
the check stars were GSC 8984:1473 and GSC 
8984:1623. Weather and time constraints did not 
permit the determination of transformations to 
the standard system. There were 1973 frames 
taken in B and 1962 in V. However, wispy clouds 
often caused scattered measurements. In order to 
remove errant measurements in an objective way, 
the raw comparison magnitudes were plotted ver- 
sus time. Points adjacent to rapidly changing val- 
ues (greathcr than 0.2 magnitudes per sampling 
interval of 30 seconds) were rejected. In this way, 

Available at http://www.axres.com 



the data set was reduced to 1148 points in B and 
1493 points in V, and the light curves were cleaned 
up considerably. 

4. Data Analysis 

We performed a simultaneous analysis of the 

AG75 photometry and radial velocities, the Hip- 
parcos photometry, and our new photometry and 
velocities with the 2003 version of the Wilson- 
Dcvinncy (WD) code. This version of WD has the 
capability of modeling the radiation of the stars by 
Kurucz (1993) atmospheres and we used this fea- 
ture in our solutions. 

In early exploratory fits with the light curve 
program of WD, it became clear that the AG75 
photometry, the Hipparcos photometry, and our 
new photometry were not consistent in terms of 
the depths of the eclipses. The AG75 and Hippar- 
cos photometry matched well, but our new pho- 
tometry had slightly deeper eclipses, requiring an 
inclination about one degree higher. Rather than 
indicating some real change in the inclination of 
the system, the difl:erences result from the fact 
that the older photometry suffers from third light 
contamination. The Digitized Sky Survey image of 
the fleld shows a companion about 15" southeast 
of TU Mus. Our aperture photometry was done 
with a 10" aperture and therefore did not include 
the companion. The effective aperture for the Hip- 
parcos observations was about 30" (ESA 1997), 
thus including the companion. The AG75 data 
also included the companion since their aperture 
was 30" . The Guide Star Catalog lists the bright- 
ness of the companion as 13™0±0'"4. In our light 
curve solutions below, we allowed for third light 
in all of the photometry but the values for the 
BV curves were always very tiny, thus we conclude 
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that third hght does not affect our BV data. The 
solution did find third light in the AG75 and Hip- 
parcos data and the results are consistent with the 
GSC brightness of the companion. Figure 5 shows 
the DSS image of the area around TU Mus. 

Previously, WD had been modified to allow for 
the use of either binary phase or time as the in- 
dependent variable (Wilson & Terrell 1998). In 
the latter case, one can solve for the orbital pe- 
riod (P), its first time derivative (P), and the 
reference epoch {HJDo, usually the time of pri- 
mary minimum). In order to investigate the be- 
havior of the period of TU Mus, we used time as 
the independent variable for our simultaneous so- 
lutions, adjusting the period and reference epoch. 
Attempts to adjust P always resulted in values 
statistically indistinguishable from zero, indicat- 
ing that the period of TU Mus has remained con- 
stant over the thirty year span of the observations 
being analyzed. Our estimate of the ephemeris is 

HJDmin = 2448500.3080(1) + 1.38728653(2) x E 

(2) 

for primary eclipse where the quantities in paren- 
theses are the Icr uncertainties in the last digits of 
the parameters. 

Other parameters adjusted in the simultaneous 
solution were the semi-major axis of the relative 
orbit (a), the binary center of mass radial velocity 
(1^), orbital inclination (i), secondary mean effec- 
tive temperature (T2), common envelope surface 
potential (f2), mass ratio (g), and the bandpass- 
specific luminosity of the primary (ii). Other 
parameters, such as the bolometric albedos and 
gravity brightening exponents, were held fixed at 
their expected theoretical values. The square root 
limb darkening law was used with coefficients from 
Van Hamme (1993). The mean effective temper- 
ature of the primary was set to 35,000 K based on 
the 07.5 spectral typo. Data set weights were de- 
termined by the scatter of the observations. WD 
gives output at each step from which the standard 
deviations of the data sets, tr, can be calculated. It 
is sufficient to determine reasonably good starting 
values for the standard deviations by examining 
plots of the observations and then adjust them 
as the solution progresses. Within a given light 
curve or radial velocity curve, WD can also apply 
individual weights. For the radial velocities, we 
set the base curve weights for the AG75 velocities 



and then gave our new velocities higher individ- 
ual weights based on their estimated errors. In 
order to eliminate potential local minimum prob- 
lems, we began the solutions from several starting 
points and found that the same minimum was al- 
ways recovered. 

Table 4 shows the results of the simultaneous 
solution. TU Mus is just barely in an overcontact 
configuration, having a surface potential of 3.137 
± 0.002 while the inner critical potential is 3.156 
for the estimated mass ratio of 0.651. Using the 
degree of contact index defined by Wilson & Rafert 
(1981), where point contact has a value of 1.0, we 
find a value of 1.006 for TU Mus. Figure 6 shows 
the fit to the AG75 radial velocities and the ones 
we obtained. The fit to the Hipparcos photometry 
is shown in Figure 7 and since all of the various ob- 
servations were phased with the ephemeris given in 
Equation 2, it shows that the period has remained 
constant to the precision with which we can mea- 
sure it. Figure 8 shows the fit to the AG75 u obser- 
vations and exhibits a small but noticeable asym- 
metry in the maxima, a phenomenon frequently 
seen in the light curves of late-type overcontact 
(W UMa) systems. As can be seen in Figure 9, 
the asymmetry is slightly smaller in the y data, 
indicating that the source of the asymmetry is a 
high temperature phenomenon, perhaps a superlu- 
minous area as suggested in CE Leo by Samec, et 
al. (1993) and YZ Phe by Samec & Terrell (1995). 
Figure 10 shows the fits to our B and V data. An 
asymmetry is present in these newer data but in 
the opposite sense from the AG75 data, indicat- 
ing that the source of the asymmetry is variable 
in time. 

5. Age and Evolutionary Status 

Very little is known about the structure and 
evolution of overcontact binaries with radiative en- 
velopes. To our knowledge, no models of this type 
have been published. Although obviously limited, 

a comparison with single-star models can be in- 
structive in estimating the age and evolutionary 
status of TU Mus. Using a stellar evolution code 
(Han, et al. (1994), Eggleton (1973), Eggleton 
(1972), Eggleton (1971)) kindly supplied to us 
by P. Eggleton, we have constructed evolutionary 
models of 23.5AfQ and 15.3Af0. The model for 
the primary reaches the observed radius of the TU 
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Table 4 
Parameters of TU Mus 



Parameter 


Value"- 


a 


17.7±O.2fl0 




—4 ± 4 km sec~^ 


i 


77°8 ± 0°1 


Ti 


35,000 k 


T2 


31,366 ± 16 K 


fii 


3.137 ±0.002 


g 


0.651 ± 0.001 


HJDa 


2448500.3066 ± 0.0001 


P 


1.38728653 ± 0.00000002 days 


Li/(Li + L2)hp 


0.648 ± 0.006 


Li/{Li + L2)u 


0.659 ± 0.003 


+i2)„ 


0.648 ± 0.003 


Li/(Li + L2)b 


0.646 ± 0.003 


+ L2)» 


0.646 ± 0.003 


Li/(Li + Z,2)b 


0.648 ± 0.003 


Li/{Li + L2)v 


0.646 ± 0.003 


h{Hp) 


0.030 ± 0.007 


h(u) 


0.023 ± 0.004 


hiv) 


0.023 ± 0.004 


h(b) 


0.025 ± 0.004 


hiv) 


0.026 ± 0.004 


Ri 


7.48 ± 0.08 Rq 


R2 


6.15 ± 0.07 Rq 


Ml 


23.5 ± 0.8 M0 


M2 


15.3 ± 0.4 Mq 


log Li/Lq 


4.8 ±0.2 ■= 


log 1/2/1-0 


4.5 ±0.2 



^Quoted errors are l-tj errors. 

^Is values are in units of total system light at phase 
0.25. 

"^The luminosity errors are estimates based on the 
uncertainty of a few thousand K in the effective tem- 
perature of the primary. The error contributions to 
the luminosities due to the errors in the radii are more 
than an order of magnitude smaller. 



8 




Fig. 6. — The fit to the radial velocities from this 
paper (solid symbols) and those of AG75 (open 
symbols). 
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-0 2 0.0 0.2 0.4 

Phase 



Fig. 8. — The fit to the AG75 u observations using 
the parameters of Table 4. 




-0 2 0.0 0.2 0.4 

Phase 



Fig. 9. The fit to the AG 75 y observations using 
the parameters of Table 4. 
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1.1 




0.4 -I , 1 , 1 , 1 

-0 2 0.0 0.2 0.4 

Phase 



10. The fit to our new BV observations 
g the parameters of Table 4. 



Mus primary at about 2.6 Myr and the model for 
the secondary reaches the observed radius at about 
6.2 Myr. Because of the luminosity exchange from 
the primary to the secondary, the radius of the 
primary is probably less than it would be with- 
out the luminosity transfer. Similarly, the sec- 
ondary is probably larger than it otherwise would 
be. Thus the true age of TU Mus probably lies 
between these two values. Our ignorance of the 
structure and evolution of these stars is so great, 
however, that one cannot discount the possibility 
that the current configuration is a result of a much 
more complicated evolutionary history. 

Two results from this work, the overcontact 
configuration and the lack of any period change, 
have interesting implications for the structure of 
early- type overcontact systems. In the late- type 
overcontact (W UMa) systems, the luminosity 
transfer drives mass transfer from the less mas- 
sive to the more massive star, resulting in period 
changes. But if TU Mus truly is in an overcon- 
tact configuration, the same level of mass transfer 
cannot be taking place because of the constant 
period. What is clearly needed is more work in 
this area with modern codes like Djehuty (Eggle- 
ton, et al. 2003). More sophisticated modeling 
by 3D stellar evolution codes will certainly give us 
a better understanding of the structure and evo- 
lution of early-type overcontact systems and the 
well-determined properties of TU Mus makes it a 
good target of such modeling. 

6. Conclusions 

Our analysis of new and existing data on TU 
Mus shows it to be an overcontact binary consist- 
ing of a 23.5M0 primary and a 15.3M0 secondary 
with a small degree of contact. Our results are in 
very good agreement with the results of AG75 and 
contradict the smaller masses found by S95. The 
source of the discrepancy between the UV and op- 
tical velocities is not clear. One idea we explored 
was the effect of the change in the center of light 
due to the reflection effect between the UV and 
optical, but our experiments indicated that this ef- 
fect is much too small to explain the large discrep- 
ancy. Given our agreement with the AG75 veloci- 
ties, one is left with two possible conclusions: that 
the UV velocities suffer from some unrecognized 
systematic error or there is some fundamental ef- 
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feet that eauses a difFerenee between UV and op- 
tical velocities for stars of this type. We note that 
Stickland and colleagues have seen a similar dis- 
crepancy for LY Aur (Stickland, ct al. 1994). Our 
current understanding of the structure and evolu- 
tion of early-type overcontact binaries is quite lim- 
ited and we hope that our results for the observed 
properties of TU Mus may serve as a good test of 
future 3D models of these stars. 
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era. RHN would also like to thank specifically the 
two technicians at the observatory, Alan Gilmore 
and Pam Kilmartin, for their superb help, and 
the members of the Physics and Astronomy De- 
partment at the University of Canterbury for their 
warm hospitality. This research has made use of 
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